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Обосновывается применение фазочастотных характеристик отраженных волн в качестве информативных параметров для реше- 
ния задач прогнозирования геологического разреза. Приводятся результаты исследования эффективности разработанных фазо- 
частотных методов анализа сейсмических волновых полей. Показана высокая разрешающая способность и помехоустойчивость 
этих методов. Выданный на основе анализа полевых сейсмических материалов прогноз нефтегазоносности подтверждается ре- 
зультатами бурения. 


В настоящее время при анализе волновых сейс- 
мических полей широко применяют динамические 
параметры волн, связанные с амплитудой и энер- 
гией отражений. В значительно меньшей степени 
используются фазовые характеристики (ФЧХ) 
отраженных волн. Между тем, в фазу сейсмических 
сигналов, а вернее, в сложный закон изменения их 
фазовых спектров, заложена важная информация о 
местоположении отражающих границ, типе ско- 
ростного разреза, поглощающих и дисперсионных 
свойствах слоистых сред [1—5]. На этой основе мо- 


гут быть синтезированы фазочастотные алгоритмы 
обработки сейсмических данных, позволяющие в 
условиях априорной неопределенности относи- 
тельно формы исследуемых волн обнаруживать и 
разрешать сигналы на фоне интенсивных помех и 
получать надежные оценки их параметров. Ис- 
пользование фазочастотных характеристик отра- 
женных волн оказывается также полезным при ре- 
шении задач прогноза геологического разреза 
[1—5], в том числе прогнозе нефтегазоносности ос- 
адочных толщ. В статье кратко обобщаются резуль- 
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таты разработки и исследования ряда фазочастот- 
ных методов анализа сложных волновых полей и 
даются отдельные примеры их применения для ре- 
шения задач прогноза геологического разреза. 

Прослеживание сейсмических волн является 
одной из центральных задач структурной сейсмо- 
разведки. Эта процедура имеет важное значение и 
для решения задач прогноза геологического разре- 
за, в том числе прогноза залежей углеводородов. 
При решении этих задач не только требуется оце- 
нивать время прихода целевых волн, но и осущест- 
влять детальное расчленение слоистых неоднород- 
ных осадочных толщ, определять форму выделя- 
емых геологических объектов. В [3] предложен оп- 
тимальный метод прослеживания фиксированных 
волн на фоне гауссовых помех, который реализует- 
ся как процедура поиска положения максимума 
функции вида: 

С(т) = ^у(/ к )со$(А(р(/ к )-2к/ к т). (1) 

к=\ 

В выражении (1) весовая функция у(/ к ) опреде- 
ляет отношение сигнала к шуму, а А (р(/ к ) - разность 
фазовых спектров сигнала и смеси «сигнал плюс 
шум» на частоте у*, т - сдвиг окна анализа, п - чи- 
сло частот, к=1,...,п. 

При практическом применении оптимального 
фазочастотного метода прослеживания (ФЧП) воз- 
никают определенные трудности, связанные, в част- 
ности, с оценкой величины уі/ к ) в выбранном диапа- 
зоне частот. Поэтому в [5] были предложены фазоча- 
стотные алгоритмы, названные алгоритмами с рав- 
новесной и неравновесной обработкой. Эти алгорит- 
мы получены из оптимального метода путем замены 
в (1) весовой функции у(/' к ) на другие, специально по- 
добранные функции. В общем случае, функция ка- 
чества Цт) для таких алгоритмов фазочастотного 
прослеживания может быть выражена в виде: 

Цт) = ^Щ/ к )со8 (<р(/ к ) —2л/ к т), 

к = 1 

где Щ/ к ) - принимаемая весовая функция, (р(/ к ) - 
мгновенный фазовый спектр выделяемого участка 
сейсмотрассы. Временное положение волн оцени- 
вается по положению максимума функции Цт), 
которая формируется при перемещении окна ана- 
лиза г вдоль сейсмической трассы, рис. 1. Для рав- 
новесного алгоритма весовая функция Щ/ к ) при- 
нимается равной единице во всей используемой 
полосе частот; для алгоритмов с неравновесной об- 
работкой - задается треугольной, синусоидальной 
или экспоненциальной формы. 

Как показали статистические исследования 
[3,4], при достаточно широкой полосе анализиру- 
емых частот и отношениях сигнала к шуму больше 
единицы данные алгоритмы незначительно прои- 
грывают в точности оценки местоположения сиг- 
нала оптимальному фазочастотному алгоритму, 
рис. 2. 



Рис. 1. Схема фазочастотного прослеживания сейсмических 
волн; т ~ текущее положение центра окна анализа, 
Т ~ длительность окна анализа, ( ~ время 



Р 

Рис. 2. Зависимость среднеквадратичного отклонения о 
оценки положения сейсмосигнала от отношения сиг- 
нала к шуму р. Алгоритмы: 1) с равновесной обработ- 
кой, 2) оптимальный фазочастотный 
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Рис. 3. Оценка разрешающей способности алгоритма ФЧП с 
неравновесной обработкой: а) модель сейсмограм- 
мы; б) разрешающая способность алгоритма: АН ~ 
оценка разрешения по вертикали; А - длина волны 


При решении задач детального расчленения 
тонкослоистых продуктивных толщ особо важное 
значение приобретает обеспечение повышенного 
разрешения выделяемых сейсмических сигналов. 
На рис. 3 в качестве примера показаны результаты 
исследования разрешающей способности алгорит- 
мов ФЧП с неравновесной обработкой, получен- 
ные на модели поля, включающей две одинаковые 
по форме волны. Как видно из рис. 3, данный алго- 
ритм обеспечивает разрешение по вертикали по- 
рядка четверти длины волны. Предложенные алго- 
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ритмы ФЧП используют для прослеживания фик- 
сированных волн, а также для построения объект- 
но-ориентированных разрезов, которые отличают- 
ся повышенной детальностью и разрешением. 

При решении задач прогнозирования геологи- 
ческого разреза важным вопросом является изуче- 
ние взаимосвязи динамических характеристик 
отраженных волн с петрофизическими параметра- 
ми геологических сред. Один из широко распро- 
страненных подходов к построению моделей 
слоистых поглощающих сред связан с представле- 
нием среды в виде линейной системы, вносящей 
определенные изменения в проходящие через нее 
колебания [1]. В рамках такого подхода было про- 
ведено исследование влияния акустической же- 
сткости и поглощающих свойств сред на фазовые 
спектры отраженных волн [6]. Исследования про- 
водились для моделей поглощающих сред, состоя- 
щих из ряда однородных слоев с плоскопараллель- 
ными границами при нормальном падении пло- 
ских волн. Спектральная характеристика системы, 

Модель 1 

Аргиллиты 

Газонасыщенный песчаник 
Аргиллиты 



описывающая свойства среды, может быть пред- 
ставлена как: 

н^Л=к м М)+ 



здесь М - количество слоев; Н т (/) — частотная ха- 
рактеристика /77-го слоя; Кіі/) - комплексный ко- 
эффициент отражения на кровле /'-го слоя; г тт _ х (]) 
— комплексный коэффициент преломления при 
прохождении волны из т-го слоя в /77-1-й. 

В качестве примера на рис. 4 показаны полу- 
ченные в результате исследования 4-х моделей по- 
глощающих сред обобщенные фазовые спектры 
волн, отраженных от подошвы 2-го слоя. Эти моде- 
ли имитировали наличие слоя песчаника, располо- 
женного во вмещающей толще аргиллитов. В пер- 
вой модели песчаник предполагался насыщенным 
газом, во втором - нефтью, в третьей - водой. В че- 
твертой модели имитировалось наличие карбони- 
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Рис. 4. Обобщенные фазовые спектры отраженных волн для принятых моделей поглощающих сред: 7 0 - базовая; і ~ текущая 
частота 
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зированного песчаника. Из графиков, предста- 
вленных на рис. 4, видно, что насыщение породы- 
коллектора газом или нефтью приводит к резкому 
изменению фазовых спектров волны, причем в ос- 
новном диапазоне частот их значения сосредоточе- 
ны в области, прилегающей к -я/2. В случае водо- 
насыщенной песчаной толщи фазовые спектры 
слабо зависят от частоты, и их значения попадают в 
область — (2,5...3) рад. При отсутствии коллектора 
(модель 4) значения фазового спектра волны, отра- 
женной от подошвы 2-го слоя, близки -я рад. 

Проведенные на ряде моделей геологических 
сред исследования показали, что фазовые спектры 
отраженных волн значимо зависят от изменения 
их акустических свойств. Это является важной 
предпосылкой для использования этих спектров в 


качестве диагностических признаков при прогнозе 
геологического разреза осадочных толщ. 

На основании данных исследований были раз- 
работаны алгоритмы фазовременного (ФАН) и фа- 
зоэнергетического (ФЭА) анализов. ФАН-анализ 
используется при прогнозе типов геологического 
разреза в около- и межскважинном пространстве 
[6]. Для каждой анализируемой сейсмотрассы в 
области исследуемой толщи формируется набор 
функций качества ФЧП общего вида: 

ь. (т) = Ж. (/ к ) С 08 (<р(/ к )-2к/ к т), (2) 

к = 0 

где ] - номер фильтра (/=1 число использу- 
емых фильтров. Совокупность функций вида (2) 
образует ФАН -колонку (ФАН -образ) для данной 
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Рис. 5. Фрагменты временных разрезов: а) ОГТ и 6) ФЧП одного из профилей Мыльджинского газоконденсатного месторож- 
дения. Отмечены положения меловых (Б & ~Б 20 ) и верхнеюрских (Юф 2 ~Ю 2 ) пластов, Вд ~ баженовская свита 
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сейсмотрассы, а их совокупность для набора сейс- 
мотрасс вдоль профиля наблюдения образует 
ФАН-разрез. 

Для решения этой задачи привлекают также 
данные геофизического исследования скважин и 
имеющиеся материалы по пробуренным ранее на 
исследуемой площади скважинам. На основе этих 
данных при проведении ФАН -анализа для пикетов 
наблюдения, соответствующих положению сква- 
жин, получают ФАН-образы, которые далее ис- 
пользуют в качестве эталонных при прогнозе типов 
геологического разреза в около- и межскважинном 
пространстве. Для проведения прогноза привлека- 
ют алгоритмы классификации, в основу которых 
положена та или иная процедура распознавания 
образов. 

При проведении ФЭА формируют функцию, 
пропорциональную энергии функций качества (2), 
получаемых при фазовременном анализе. Анома- 
лии в поглощении и дисперсии скорости сейсми- 
ческих волн приводят к появлениям аномалий в 
поведении этой функции, что может служить прог- 
нозным признаком для выделения нефтегазоперс- 
пективных зон. По данным ФАН и ФЭА-анализов 
строят карты распространения типов разрезов и ра- 
спределения аномальных зон для изучаемой пло- 
щади. По результатам комплексной интерпретации 


всех материалов фазочастотной обработки форми- 
руют предложение по размещению промысловых и 
разведочных скважин. 

На базе вышеизложенных фазочастотных мето- 
дов были разработаны программно-алгоритмиче- 
ские комплексы «Геосейф» и «Геосейф-ВСП» [5,6]. 
Комплекс «Геосейф» ориентирован на обработку 
временных разрезов, получаемых при проведении 
наземных исследований методом общей глубинной 
точки (ОГТ). Комплекс «Геосейф-ВСП» использу- 
ется для обработки данных вертикального сейсми- 
ческого профилирования скважин [6]. На этих 
программных комплексах проведена обработка 
большого объема материалов сейсмических иссле- 
дований ОГТ и вертикального сейсмического про- 
филирования, полученных в условиях Томской и 
Тюменской областей. На рис. 5-7 приведены от- 
дельные результаты обработки материалов назем- 
ной сейсморазведки с помощью разработанных 
фазочастотных алгоритмов. 

Так, на рис. 5 представлены фрагменты разреза 
ОГТ и полученного из него объектно-ориентиро- 
ванного разреза (разреза ФЧП) для одного из про- 
филей Мыльджинского газоконденсатного место- 
рождения Томской области. Как видно из рисунка, 
на разрезе ФЧП, в отличие от ОГТ, удается выде- 
лить все основные отражающие горизонты нижне- 



III - тип ГѴ - тип 

Рис. 6. Эталонные ФАН-образы, выделенные для различных типов геологического разреза. Сг ~ георгиевская, Тт ~ тюменская свиты 
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го мела и верхней юры; уверенно прослеживаются 
даже очень слабые по интенсивности отражения; 
более отчетливо проявляются различные неодно- 
родности исследуемого геологического разреза. 



Рис. 7. Проявление нефтяной залежи на фазоэнергетиче- 
ском разрезе (Болтная площадь, Томская обл.) 
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